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*\ Abstract. Superfluid properties of electron- hole pairs in a quantum Hall four-layer 

system are investigated. The system is considered as a solid state realization of a 
two-component superfluid Bose gas with dipole-dipole interaction. One superfluid 
component is formed in the top bilayer and the other component - in the bottom one. 
We obtain the dispersion equation for the collective mode spectrum and compute the 
critical parameters (the critical interlayer distance and the critical currents) versus 
the filling factor. We find that the critical currents of the components depend on 
O , each other. The maximum critical current of a given component can be reached if 

the current of the other component is equal to zero. The non-dissipative drag effect 
between the components is studied. It is shown that in the system considered the drag 
V/ factor is very large. Under appropriate conditions it can be about 10 per sent, that is 

\ jr} • at least in three order larder than one predicted for two-component atomic Bose gases. 
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1. Introduction 

Among the objects that demonstrate Bose-Einstein condensation or superconductivity 
considerable attention is given to two-component systems. Particularly, beginning from 
the paper by Andreev and Bashkin [1] the possibility of a non-dissipative drag between 
superfluid (superconducting) components moving with different velocities was discussed 
[21 [3j HI El [6] (see also the review [7]). The Andreev-Bashkin effect was also considered 
in astrophysics in the context of superfluid models of neutron stars [SI E]. A related 
problem - the critical velocities in two-component superfluid systems was studied in a 
recent paper [10] . It was shown that critical velocities are essentially different in the 
case when two components move with the same velocities and in the case when one 
of the components does not move. In the latter case the critical velocity of a moving 
component can be much higher. 

Although two-component superfluid atomic Bose gases have been realized in 
laboratories [H] , there are certain problems in experimental observation of the effects 
caused by relative motion of the components. On the one hand, it is not so simple 
to create a relative flow of superfluid components in a mixture of Bose gases. On the 
other hand, a spatial separation of components takes place in two-component mixtures 
confined in a trap. In the gases with point interaction the spatial separation results 
in a disappearance of the effects caused by inter-specie interaction. To overcome these 
difficulties one can turn to systems where components can be kept spatially separated 
in a controllable way (that gives a possibility to provide a flow of the components with 
different velocities). The interaction between the components in such systems should 
contain a long-range part. For example, one can deal with a Bose gas with dipole-dipole 
interaction confined in a double-layer trap [5]. 

In this paper we study a solid-state system where a two-component Bose gas with 
the dipole-dipole interaction can be realized. We consider a multilayer electron system 
where electrons from one layer couple with holes from the adjacent layer. Since such 
Bose particles (electron-hole pairs) have a small mass they may demonstrate superfluid 
behavior at rather high temperatures (much higher than required for the Bose-Einstein 
condensation of alkali metal vapors). Electron-hole pairs in such systems have large 
dipole momentum and the dipole-dipole interaction determines, in the main part, 
collective properties of a gas of such pairs. 

To be more specific, we consider a four-layer electron system in a strong 
perpendicular to the layers magnetic field (a multilayer quantum Hall system). For 
bilayer quantum Hall systems with total filling factor equal to unity the theory predicts 
[T2] the existence of a superfluid condensate of indirect excitons in the systems. An 
indirect exciton corresponds to a bound state of an electron belonging to one layer 
and a hole (an empty state in the lowest Landau level) belonging to the other layer. 
This prediction was partially confirmed in experiments JT3J HH [15]. Bose-Einstein 
condensation of metastable (optically generated) indirect excitons in zero magnetic field 
was also observed [T6l ITT] . As was shown in [18], multi-component excitonic superfluid 
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condensates can be realized in multilayer quantum Hall systems with even number of 
the layers and the average filling factor per layer equal to one-half. According to [18J, 
electron- hole pairs emerge in separate bilayers, i.e., a given excitonic component belongs 
to a given double-layer complex. Here we consider a four-layer system with the filling 
factors of the layers v\ = u^, ^2 = 1 — vt, ^3 = v&, ^4 = 1 — v&- In such a system one 
specie is formed by coupled electrons and holes from the layers 1 and 2, while the other 
specie - by coupled carriers from the layers 3 and 4 (Figured]). 

In the system considered a flow of electron-hole pairs is equivalent to two oppositely 
directed electrical currents in the adjacent layers. Therefore, a superfluid state of such 
pairs can be considered as a specific superconducting state. Counterflow supercurrents 
can carry an electrical current from the source situated at one end of the system to the 
load situated at the opposite end (if the interlayer tunnelling is negligible small, the 
dissipation is negligible, as well [19]). Experimentally one can provide separate contacts 
to each layer. It allows to control and measure the supercurrents in each bilayer complex. 
Therefore, it is more appropriate to formulate the problems not in terms of superfluid 
velocities, but in terms of supercurrents. 

In this paper we address two problems. First, we investigate critical currents in 
a four-layer quantum Hall system. We show that the behavior of critical currents 
is qualitative the same as of critical velocities in two-component Bose gases (with 
the important advantage that the effect in multilayers can be registered by electrical 
measurements). Then, we consider the non-dissipative drag between the components 
and compute the drag factor. We predict very large drag factor for quantum Hall 
multilayers: under appropriate conditions it can reach 10 per sent (in three order higher 
than the most optimistic estimates for atomic Bose gases). 

Our study is based on the analysis of collective mode spectra. We follow the 
approach proposed in [20] for the study of bilayer systems with zero imbalance of filling 
factors (v\ = V2 — 1/2). In Secj2]we extend the approach [20] for the case of an arbitrary 
imbalance. In Sec]3]we obtain the spectra of excitations for the four-layer system, find 
the critical interlayer distance versus the filling factor, and obtain the relation between 
the critical currents. The non-dissipative drag effect is considered in SecJU 

2. The approach 

Let us begin with the discussion of a mechanism that determines critical supercurrents in 
quantum Hall bilayers (multilayers). In bulk superconductors the restriction on the value 
of the supercurrent emerges from the requirement that the magnetic field produced by 
electrical currents should be lower than the thermodynamic critical magnetic field. For 
thin films the critical magnetic field is higher than for bulk superconductors, and at small 
thickness w of the film it increases by the law H c ex 1/w. Due to almost two-dimensional 
character of conducting layers in quantum Hall systems the critical magnetic field should 
be high. The critical current is determined by an essentially different mechanism. 
This mechanism is just the generalization of the Landau mechanism of destruction of 
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superfluidity (governed by the Landau criterion). The critical current can be found from 
the requirement that in a superconducting state the energies of all collective excitations 
being real valued and positive. Going ahead we note that electrical currents close to 
critical ones produce magnetic fields much smaller than the terrestrial magnetic field. 
And, indeed, the magnetic mechanism of destruction of superconductivity is irrelevant 
for the quantum Hall bilayers (multilayers). 

For more transparency, we will describe the approach with reference to a bilayer 
system. Let a double-layer electron system is situated in a perpendicular to the 
layers magnetic field B. The electron density p = p\ + p 2 satisfies the condition 
z/tot = 2nl 2 p = 1, where I = ^hc/eB is the magnetic length, i.e. the filling factor 
of the layer 1 is V\ — v {y < 1), and the filling factor of the layer 2 is z/ 2 = 1 — v. 
The value v is the parameter of the problem. Since z/ l5 z/ 2 < 1 all the carriers belong to 
the lowest Landau level (under assumption that the Coulomb energy is small compared 
to the energy gap between the Landau levels). We take the Hamiltonian in the lowest 
Landau level approximation: 

H = 2S 5Z ^2 v ^'(q) \ Pn(q)/v(-q) - <5„,,„'p„(0)exp ( — — J I , (1) 

n,n'=l q ^ \ / ) 

where S = L x L y is the area of the layer, V ntn /(q) = 27ce 2 e~ dq \ n ~ n /eg is the Fourier- 
component of the Coulomb potential, d is the distance between the layers, e is the 
dielectric constant, and 

Pn(q) = J2 c n( X + ^niX - «£) exp(ig x X - t?) (2) 

X 

is the Fourier- component of the electron density in the n-th layer. In (jSJ) c+(X) and 
c n (X) are the creation and annihilation operators for electrons in the n-th layer in the 
state described by the wave function ipx( T ) = exp (—iXy/l 2 — (x — X) 2 /2l 2 ). We imply 
that the interlayer tunnelling amplitude t is much smaller than the Coulomb energy 
E c = e 2 j el and neglect the tunnelling in the Hamiltonian ([1]). We will discuss the 
validity of such an approximation in more details at the end of this section. Here we 
just mention that in the bilayer systems used in experiments the tunneling amplitude 
t ~ 50 pK [13] that is in 6 order smaller than the Coulomb energy. 

The state with electron-hole pairing can be described by a BCS-like many-body 
wave function |\l/) = Ylx [ u x + ^x^i~P0 c :tP0] \vac), where hf is the creation operator 
of the hole in the layer 1 and \vac) is the vacuum state defined as a state with 
completely filled layer 1 and empty layer 2. The u — v coefficients satisfy the condition 
l^x) 2 + l^x) 2 = 1- This function can be presented in another equivalent form 



w=n 



X 



cos %c+(X) + e lVx sin %c+(X) 



|0>, (3) 



where 9x = &{X) and (px = (f{X) are arbitrary functions. The quantity 6 can be 
connected with the local filling factors ^x,i(2) = (1± cos9 x )/2. One can see that ipx 
is the phase of the order parameter Ax = (^| cf x C2x\^f) = e lipx a/^x(1 — vx) which 
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corresponds to the electron-hole pairing. At 9 = const and ip = the function ([3]) in 
the coordinate representation coincides with the famous (1,1,1) Halperin wave function 
(see, for instance, [2T]). 

In the state ([3]) the energy of the system reads as 

E = ^- J2 { i H ( X ~ X ') ~ F s( X - X ')] cos# x cos9 x > 

y X,X' 

-F D (X-X')sin6'xsin6 l X /cos(v9x -<Px>)\, (4) 

where the quantities 



e 



00 i _ e -kM 



Fs{ x ) = e - e -^ r* d « ^ 



■J/- 



2^ i-oc y/q 2 + X 2 // 4 

F DW = ^e-S f «* e^* (5 ) 



2e ./_«, y/q* + X 2 /l* 

describe the direct Coulomb interaction, the exchange interaction in a given layer, and 
the exchange interaction between the layers, respectively 

We consider excitations above a homogeneous state with a stationary superflow 
of electron-hole pairs along the x direction. Such a state corresponds to 9x = So 
independent of X and the phase tpx — QX linear in X. The energy of the homogeneous 
state is found from Eq.(j4j) and reads as 

£ (0) = 4^2 ( [W(0) " ^s(O)] cos 2 9 - r D (Q) sin 2 0„) . (6) 

Here the functions displayed calligraphically indicate the Fourier-transforms defined as 

/oo 
dXexp(— iqX)A(X). The explicit expressions for the quantities 
-oo 

in (El) are 



e 2 _^l-e- d|<?l 

,2 roo 



HW = ie -4-i^^, (7) 



„2 roo 2 2 

^b(g) = tt dk e-—J (kql 2 ) e- kd . (9) 

Fluctuations over the stationary state can be parametrized as m z (X) = cos#x — 
cos#o an d fix — ¥x — QX. The energy of fluctuations in quadratic approximation has 
the form 

S v r i 

E& = -^ Yl [™^(-?)^(#^) + 1m z {-q)K Zip {q)(p{q) + <p(-q)fC vv> (q)<p(q) . (10) 

g 
In ( TTUj) the Fourier-components of the fields m z (X) and (px are defined as 

m^) = ^J]m 2 (X)e-^, 0(g) = ^- 5>(X)e"^. (11) 



Critical currents and non-dissipative drag in quantum Hall multilayers 6 
In Eq. ( 1101) the components of the matrix K, read as 

K„{q) = H(q) - MQ) + MQ) + (MQ) ~ • FD(g + g) ^ D(9 ~ g) ) cot 2 9 , (12) 

rC zip = -icos0 o , (13) 



fc V tp(q) = sin 2 8 Q 



r M) - *>(« + « + *>(« -Q)l . (14) 

One can note that Eq. (fl2l) diverges at 9q = and 8q = tt and the approximation 
([TO]) violates. But such 9$ correspond to filing factors v\ = and z/2 = 1 (or vice versa). 
At such filling factor the density of electron-hole pairs is equal to zero and one cannot 
speak about the spectrum of collective excitation in the gas of the pairs. The cases of 
6*0 close to or tt corresponds to low density of the pairs and the approximation (TTOj) is 
valid under condition that the density fluctuation are small as compared to the ground 
state density. The latter condition is equivalent to |m 2 (g)| -C sin#o- 

To quantize the energy ( TTOl) one notes that rh z and (p are the conjugated quantities 
and the commutator of the operators that correspond to these variables is equal to 

2vr/ 2 
[m z {q)Mq')\ = -2i^r<W ( 15 ) 

The operators m z (q) and 0(q) can be expressed in terms of Bose creation and 
annihilation operators in a common way 

m z (q) = A(b q + bt q ), 0(q) = iB{b q - b± q ), (16) 

where the amplitudes satisfy the condition AB = 2ti1 2 /S. Replacing the variables 
m z (q) and (p(q) in (TTOl) with the operators (TTBT) and requiring vanishing of the terms 
containing two creation (two annihilation) operators one finds the explicit expressions 
for the amplitudes 



As a result one obtains the Hamiltonian for the collective excitations 



n* = y £E(q)(b + (q)b(q) + ^ 



where 



E{q) = 2 JJC w (q)IC zz (q) + K z Jq) (19) 



(with JC zip (q) = iK, zip (q)) is the spectrum of collective excitations. 

It is instructive to compare the spectrum (fT9l) with the Bogolyubov spectrum. Let 
us introduce formal notations e, v, 7 and n and present Eq. ( Tl9l) in the Bogolyubov 
form 



E(q) = v /e(e + 2 1 n) + hqv. (20) 
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In f )20p the kinetic energy is defined as 

e = 2Jd(Q) - F D (q + Q)-F D ( q -Q). (21) 

In the long wave limit (q, Q -C Z -1 ) this quantity is reduced to the standard expression 
for the kinetic energy e = h 2 q 2 /2M, where 

„ 2 £ fi 2 r/7/ d 2 \ fd 2 \ t /d\ dV' 

M = -pt [v 2 1 1 + p) exp [w) erfc (ts) - 1\ (22) 

is the magnetic mass of the pair (see, for instance, [22]) . 

In (1201) the density is defined as n = z/(l — v)/2irl 2 . In the limit v — > (or 
(1 — z/) — > 0) this quantity coincides with the density of electron-hole pairs (in this 
case one can easily mark out the pairs from the background, see Figure [2]). The factor 
v{\ — v) appears due to the electron-hole symmetry of the problem. 

The superfluid velocity in (J20|) is defined as 

^■Fp( g + Q)-^ D ( g -Q) 
nq 

At q, Q <C l^ 1 it reduces to the expression v = (HQ/M)(1 — 2v) that differs from the 
common expression for the superfluid velocity by the factor 1 — 2z/. At small v this factor 
approaches unity and the difference disappears, but at zero imbalance [y = 1/2) it is 
equal to zero and the last term in (1201) vanishes. This feature can be understood from 
the following arguments. The sign of Q (gradient of the phase) determines the direction 
of the current. If one describes the supercurrent as motion of electron-hole pairs then 
the direction of the current depends on the direction of the superfluid velocity and on 
the direction of polarization of electron-hole pairs (Figure [2]). Therefore at given Q the 
filling factors v < 1/2 and v > 1/2 corresponds to opposite directions of superfluid 
velocities. The factor (1 — 2u) changes its sign under substitution v —* 1 — v and its 
appearance in the expression for v reflects to the electron-hole symmetry of the problem. 
The interaction parameter 7 in (1201) is given by the expression 

Fv(q + Q)+Fv(<l-QY 



7 = 8vr/ 2 



H(q)-F S (q) + 



(24) 



2 

The first term in ( |24l) is caused by the direct Coulomb interaction between dipoles, two 
other terms correspond to the exchange interaction. One can see that at small q the 
direct interaction term reduces to 70 = 4ire 2 d/e. It is just the interaction parameter 
for a two-dimensional gas of classical dipoles (polarized perpendicular to the layer). We 
also note that in the limit d —>■ (in which the pairs do not interact with each other) 
the quantity ( 1241) approaches zero, and the collective excitation spectrum ( 1201) turns to 
the spectrum of free particles E(q) = e = h 2 q 2 /2M (at Q = and ql <^ 1). 

Thus, the spectrum (TT9l) is similar to the Bogolyubov spectrum but there are certain 
differences caused by the electron-hole symmetry. The answer (TT9l) at v — 1/2 coincides 
with the spectrum obtained in [20]. In the low density limit v ^C 1 (low concentration 
of the pairs) and at small q the result ( TT9l) reduces to one found in |23j on the base of 
the Gross-Pitaevskii equation. 
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According to the Landau criterium the energy (j!9p should be positive and real- 
valued for all q. At d > d c (critical d c depends on v) the spectrum becomes complex- 
valued already at Q = 0. It means that at given d and v the superfluid state cannot be 
realized at all. At d < d c the Landau condition violates at Q > Q c . The value Q c (that 
depends on v and d/l) determines the critical current. The relation between Q c and 
the critical current can be obtained as follows. At nonzero vector potential the phase 
gradient Q in the energy (jfJJ) should be replaced with the gauge invariant quantity 

Q^Q-(e/hc)(A 1 -A 2 ), (25) 

where Aj is a vector potential in the z-th layer. In a given layer the current (density 
of the current) is obtained from the equation j n = cde/dA n , where e = E^ / S is the 
energy per unit area. Taking into account (1251) . one finds ji = —ji = —(e/hS) dE^/dQ, 
or, in the explicit form, 

el . 2 d.F D (Q) 
Jl = - j2 = hA^ Sm6o ^Q- (26) 

Substituting Eq.© into (|26|) one obtains the following relation between Q and the 
current 



3\ = -32 



r°° ~ k 2 

1/(1-1/)/ dke--~U 1 (kQl)e- kd/l . (27) 

Jo 



27rel 2 H 

At small Q (Q -C /) the current is proportional to the phase gradient: j\ = (e/h)p$Q, 
where p s = h 2 n/M is the superfluid stiffness. 

One can easy check that the value of integral in (|2"T|) at any Q and d does not 
exceed 0.45. Therefore, in any case the supercurrent is less than j m w 0.018e 3 /£:/ 2 /i. At 
typical parameters e = 12.5, I = 100A one evaluates j m ~ 5 A/m. Note this value of 
counterflow currents corresponds to very small value of parallel to the layers component 
of magnetic field produced by the current (B y ps 6 • 10 -6 T). Under accounting the 
Landau criterium this quantity is even smaller. 

Now let us discuss the influence of tunnelling on the critical current. In general 
case the tunnelling reduces the critical current, moreover it may destroy the superfluid 
state completely [20]. To explain this effect it is instructive to introduce the Josephson 

length A = l\/27rp s /ty/i/(l — v) (see, for instance, [E],[20]). This quantity determines 
the size of the soliton (Josephson vortex) in the bilayer quantum Hall system. If the 
Josephson length is small A < 2n/Q c and the gradient of the phase caused by the soliton 
is large (larger than the critical gradient Q c ) the tunnelling influences significantly on 
the critical parameters. But such A correspond to quite large tunnelling amplitudes 
(t/Ec > 10~ 2 for v = 1/2 and d ps /). Thus, in a common experimental situation 
(t/Ec ~ 10~ 6 ) the tunnelling can be neglected almost for all relevant values of v and 
d/l. 
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3. The excitation spectrum and critical parameters for the four-layer 
system 

Let us turn to the consideration of the four-layer quantum Hall system where a two- 
component superfluid gas of electron- hole pairs can emerge. For definiteness, we consider 
the system with equal distances between the adjacent layers. We specify the case where 
supercurrents in both components are directed along the same axis. 

The Hamiltonian of the system has the form (CEJ) (with the summation over four 
layers). According to the result of Refs. [18] the electron-hole pairs are formed separately 
in the n = 1,2 bilayer complex and in the n = 3,4 bilayer complex. Correspondingly, 
the many-body wave function can be presented as a product of the functions ([3]): 

m = n ( cos if c+ix + sin e ~f- eiVTX ctx ) ( cos ^f- ctx + sin e -f- eiVBX 4x ) |0) - (28) 

In the state ( 1281) the energy of the system consists of three terms: 

E = E T + E B + E TB . (29) 

In (1291) Ei and E B are the bilayer energies (given by Eq. (j3j)). The cross term has the 
form 

Etb = 2T ^ Htb{X ~ X>) cos6t ^ cos ^b(A'), (30) 

y X,X' 



/oo 
dqe iqX n TB (q) with 
-oo 



p 2 n _p- rf igh 2 2,2 

n TB (q) = ^e" ^ 1 , } e'^. (31) 

Note that in the state (I2"8"j) the cross energy does not contain the exchange part. 

The stationary homogeneous state is described by four parameters: z/j = (1 + 
cos6'oj)/2 and Qi = d(fix/dX (i = T, B). The energy in this state reads as 

E<® = 4 0) + 4 0) + ^i^ TB (0) cos#ot cos^ob, (32) 

where E^' (E B ') are determined by the equation (jSJ) with 8 = #ot (#ob) and 

Q = Qt(Qb). 

The energy of fluctuations is found by the same procedure as for the bilayer system. 
The result is 

Efl = iS^{ £ [*M-?)*£(?)riM?) +2mU-g)X*(g)fc(?) 

g i=T,B 

+^(-g)/C;V(g)^(g)]+m T ,,(-g)/C™(g)m B ,,(g)}, (33) 

where diagonal in i components of the matrix /C are given by the expressions (fT2]) - (fT4"j) . 
and the non-diagonal in i component is /C™(g) = H TB (q). 
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Replacing m i)Z (q) and <fi(q) by the operators 



rh>i yZ {q) = 


l2nP 

V 5 




) 4 (fc* + K- 


0i(q) = i\ 


1*? ( 




4 



qj- 



(34) 

we obtain the Hamiltonian written in terms of Bose creation and annihilation operators 
Hr = Z [ Z (^.*(?) + ^9«<) ( b tg^)Ki + l) + 9«( b iM* + b %M-« + /i.e.)] (35) 

n i—T R \ / 



with Eo,i(q) = 2.//C* (g)£« (g) and 




-, L/ |]:, 4/ ^VyJj^VyW / Qfi N 

In ( 155]) the quantities Uj are given by Eq. (T25J (with Q = Qi and z/ = z/j). 

The Hamiltonian (1351) coincides in form with one obtained in [6] for the two- 
component superfluid Bose gas. Its diagonalization yields 

^ = EE e m ( b i« b <*« + 1) > ( 37 ) 

q a=l,2 ^ ' 

where 6+ (b a ) are the operators of creation (annihilation) of collective excitations, E a (q) 
are the excitation spectra. The dispersion equation for the spectra is analogous to one 
obtained HD1: 



[El T (q) -{E- hv T q) 2 ][El B (q) - (E - hv B q) 2 } - 4g 2 E , T (q)E 0:B (q) = 0. (38) 

One can show (see [TU]) that the energies of collective excitations are real valued and 
positive if the quantities E j (q) are real valued for all q, and the following inequalities 
are satisfied 

[El T (q) - (hv T q) 2 }[El B (q) - (hv B q) 2 } - 4g 2 q E , T (q)E , B (q) > 0, (39) 

Eo,T(q) - hv T \q\ > (40) 

(the condition (1401) can be replaced by E 0<B (q) > hv B q). 

Putting Qt = Q B = and solving Eq. (|38l) one obtains the spectra of collective 
modes at zero currents. The requirement of real valued spectra yields the critical 
interlayer distance d c . The dependence of d c on the filling factors is shown in Figure 
[3] (for z/t = z/r). The behavior of the critical interlayer distance is the same as for the 
bilayer system [24], but absolute values of d c are a little bit smaller. The minimal critical 
distance (rf c , m in ~ 1.015/) corresponds to the case of zero imbalance of filling factors. At 
d > ^c,min one can say also about critical filling factors (that decreases under increase 
oid). ' 

A state with nonzero supercurrents can be realized only d < d c (i>^, z/r) if for given 
Qi the spectra satisfy the inequalities (|39|) . (1401 . The currents are determined by the 
relation j, = —(e/hS)dE^/dQi. Since the cross term in Eq. (l32~l) does not depend on 
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Qj, the relation between Qi and ji is same (Eq. (125|) ) as for the bilayer system (if one 
neglects the fluctuating part of energy, see Sec. @J. The inequalities fl39|) . f|40|) determine 
a joint condition on Qi and the critical current of one component depends on the current 
of the other component. 

The relations between the critical currents for d/l = 0.9 and v^ = v-q are shown in 
Figure HI According to Figure HI typical absolute values of the critical currents are less 
or of order of 1 A/m. Comparing the results presented in Figure H] with ones of Ref. [TO] 
one can see that the critical currents demonstrate the behavior similar to one of critical 
velocities in two-component superfluid Bose gases. Namely, the maximum supercurrent 
of one component can be reached at zero supercurrent of the other component, while 
at equal currents their allowed values are the smallest one. Since measurements of 
electrical currents in the layers are more simple than the measurements of superfluid 
velocities in two-component mixtures, quantum Hall four-layer systems can be used for 
the observation of specific behavior of critical velocities in two- components superfluids 

ma. 



4. Non-dissipative drag between the components 

Eq. (l26l) used in the previous section for the calculation of the current does not take 
into account the energy of fluctuations. Therefore, the results obtained are valid at 
temperatures much smaller than the Coulomb energy. In two-component systems the 
fluctuations yield an additional contribution to the current even at zero temperatures. 
This contribution is caused by the energy of zero-point oscillations. It is rather 
small contribution and it can be neglected under calculations of the critical currents. 
Nevertheless, this contribution determines a new effect - a non-dissipative drag between 
the components. The value of the non-dissipative drag decreases under increase of the 
temperature, but the decrease of the drag factor is essential at temperatures larger than 
the interaction energy [6] . Here, for simplicity, we consider the case of zero temperatures 
and small phase gradients (Qil <Cl). 

Taking into account the zero-point oscillations energy 

^^ (0, + ^E^) ( 41 ) 

a=l,2 q 

and expanding it in series in Qi one obtains the following expression for the energy 

E^E + -J2 A ikQiQk, (42) 

ik 

where E is independent of Qi, and A is some symmetric real matrix. The mean-field 
energy E^ in ( 14T1) is diagonal in Qi, but the zero-point oscillation energy contains a 
non-diagonal term. Due to this a supercurrent of a given component depends on the 
phase gradients of both components 

3i = ^^^ikQk- (43) 
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The latter results in the non-dissipative drag effect. Indeed, let the current in the 
drive component (e.g. component T) is given (it is fixed by an external source) while 
the current in the drag component (B) is not fixed. The value of current in the drag 
component can be found from the requirement of minimum of the energy Eq. (1421) 
subjected to the constrain j't = const. One can see that the drag current is nonzero 
and proportional to the non-diagonal component of the matrix A: 

JB = ^jr- (44) 

We define the drag factor as the ratio of the drag current to the drive current: 
/dr — Abt/Att- As was found in [5j [6] the drag factor for the atomic Bose gases is 
rather small - the most optimistic estimates yield /dr ~ 10~ 4 . Let us compute the drag 
factor for the quantum Hall multilayers. 

The main contribution into the diagonal components of A comes from the energy 
E(°\ and the quantity Att in the leading order is evaluated as Att = f?n^jM 
irii = Ui(l — v{)/2irl 2 ). The non-diagonal component of A is caused by the zero-point 
oscillation energy. Strictly spearing, to compute this energy one should obtain the 
spectrum of collective excitations for all q (not only for q || x). But if one needs only 
non-diagonal in Qi term and Qi are assumed to be small an approximate dispersion 
equation can be used. This equation is obtained from fl38|) if one replaces the quantities 
hviq with hviq x and neglects the dependencies of E 0ti on Qi. This approximation can 
be justified under accounting that in the series for E a (q) the QtQb terms comes only 
from the product ^b^t (see [6]). 

Such an approximation yields the dispersion equation 

[E 2 -{E- hv T q x ) 2 }[El -(E- hv B q x ) 2 } - 4( 7 ') 2 e 2 n T n B = 0, (45) 



where E t = i/e[e + 2^nj\ are the spectra of excitations for decoupled one-component 
systems, e = 2 [^d(O) — ^(q)] is the kinetic energy of electron-hole pairs, 7 = 
8nl 2 [H(q) — J-"s(q) +-7"d(0)] and 7' = 4tt1 2 7{ tb (q) are the Fourier components of the 
interaction potentials, and Vi = (2/hq) (djF D (g)/dg) (2z/j — l)Qi are the superfluid 
velocities (with the factors that account the electron- hole symmetry). 

Eq. f|45|) coincides in form with one for the atomic two-component Bose gases 
P, [TO] • For obtaining Abt we present the solutions of ( )45l) as series in Vi and substitute 
them into (j4Tj) . The details of such a procedure are described in [6]. Here we present 
the final expression for the drag factor 

BT im r (V)V^b6 2 Air D \ 2 

/dr =—r 2 — l-2z/ T l-2z/ B / —————— 3 I -5— gdg, 46 

nh 2 n T J E a Ep{E a + E p ) 6 \ dq J 

where 



EaM = \l ^^ ± f El A ElY + 4(7') 2 ^ T n B (47) 

are the energies of collective excitations at Qt = Qb = 0. 
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Since the problem has one energy parameter (the Coulomb energy e 2 /el) the drag 
factor (|46p depends only on dimensionless quantities d/l, z/p and z/b. The dependence 
of the drag factor on the filling factor [y = z/p = ub) at different d/l is shown in Figure 
[5j One can see that this dependence has an extremum at small v (small density of 
the pairs). Similar feature was obtained for the atomic Bose gases with dipole-dipole 
interaction [5]. 

The specific feature is a sharp increase of the drag factor near critical v (or d, see 
Figure EI). The effect is caused by a roton-like minimum in the spectrum of the lowest 
collective mode (Figure [7]). 

Other specific features are the vanishing of the drag effect at z/b = 1/2 (z/p = 1/2), 
and the alternation of the sign of the drag factor under change of sign of the filling factor 
imbalance (z/p to 1 — i/p or z/b to 1 — z/b)- The alternation of the direction of the drag 
current can be observed if one keeps one of the filling factor constant, while the other 
filling factor tunes from z/j < 1/2 to z/j > 1/2. This feature can be understood from the 
electron- hole symmetry argumentation (see the discussion in Sec. [2]). 

But the main feature is large absolute values of the drag factor (in 10 2 -r- 10 3 times 
larger than the most optimistic figures for the atomic Bose gases). Large values of the 
effect are caused by a number of factors. As follows from the consideration [5j [6] the 
drag effect in two-dimensional systems can be larger than in three-dimensional ones, 
but for Bose gases in bilayer traps large values are not reached due to weak interspecie 
interaction. Fortunately, in quantum Hall multilayers the interaction between different 
superfluid components is of the same order as the interaction inside a given component. 
Moreover, the intra-component interaction is reduced due to the exchange interaction. 
At last, the drag effect is enhanced considerable by the presence of the roton-like 
minimum in the energy spectrum, and this minimum becomes deeper at interlayer 
distances or filling factors close to critical ones (Figure [7J . All these factors work in 
favor of the drag effect and result in a giant drag factor in comparison with atomic 
two-component Bose systems. 

5. Conclusion 

In conclusion, we have studied superfluid properties of a two-component gas of electron- 
hole pairs in a quantum Hall four-layer system. We have found that the critical 
parameters (critical interlayer distance and critical currents) for this system are slightly 
less but of the same order as in bilayers. The critical currents in the two-component gas 
of electron-hole pairs demonstrate the behavior similar to one of atomic two-component 
Bose gases. In particular, the largest value of the critical current in a given component 
can be reached if the current in the other component is equal to zero. In multilayer 
quantum Hall systems this peculiar behavior of two-components superfluids can be 
observed by electrical measurements. The non-dissipative drag effect between the 
components is predicted. The effect takes place only at nonzero imbalance of filling 
factors of each (top and bottom) bilayer, and the drag current alternates its direction 
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under change of sign of the imbalance in one of the bilayers. The drag factors is quite 
large, and the largest values can be achieved at the interlayer distances close to the 
critical ones. 
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Figure 1. The four-layer system with two-component gas of electron-hole pairs. 
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Figure 2. Relative directions of electrical currents in the layers (j) and the velocities 
of electron-hole pairs (v) at v < 1/2 (a) and v > 1/2 (b). 
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Figure 3. Critical interlayer distance for the four-layer system vs the filling factor 

{v = v T = v B ). 
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Figure 4. Critical currents (in e 3 /hel 2 units) at d/l = 0.9 and the filling factors 
vt = vb = 0.5, 0.25, 0.15, 0.07 (from the top to the bottom curve) 
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Figure 5. Drag factor vs the filling factor [y = vt = Vb) 
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Figure 6. Drag factor vs the interlayer distance (at vt = vb = 1/4). 
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Figure 7. The energies (in e 2 /e? units) of the collective modes at vt = Vb = 1/4. 
Dashed curves - d= I, solid curves - d = 1.16Z 



